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ENVIRONMENTAL FATIGUE OF AN Al-Li-Cu ALLOY: 
PART III - MODELING OF CRACK TIP HYDROGEN DAMAGE 


Robert S. Piascik and Richard P. Gangloff' 


ABSTRACT 

Environmental fatigue crack propagation rates and microscopic damage modes in Al-Li-Cu alloy 
2090 (Parts I and II) are described by a crack tip process zone model based on hydrogen 
embrittlement. Da/dN ENV equates to discontinuous crack advance over a distance, A a, determined 
by dislocation transport of dissolved hydrogen at plastic strains above a critical value; and to the 
number of load cycles, AN, required to hydrogenate process zone trap sites that Jtoaune 
according to a local hydrogen concentration-tensile stress criterion. Transgranular {100} 
cracking occurs for process zones smaller than the subgrain size, and due to lattice decohesion 
or hydride formation. Intersubgranular cracking dominates when the process zone encompasses 
one or more subgrains so that dislocation transport provides hydrogen to strong boundary 
trapping sites. Multi-sloped log da/dN-log A K behavior is produced by process zone plastic 
strain-hydrogen-microstructure interactions, and is determined by the A K dependent rates and 
proportions of each parallel cracking mode (0;) according to: 

da/dN ENV a 0{m} AK + e <ioo} AK + e sac AK * 

Absolute values of the exponents and the preexponential coefficients are not predictable , however, 
fractographic measurements of 0, coupled with fatigue crack propagation data for alloy 2090 
establish that the process zone model correctly describes fatigue crack propagation kinetics. For 

eXampk tei<m^ Al , = e |llll ( 3 xio’)AK 40 + e.^oxiow; +0(2110^* 
where each 0 is known as a function of AK. Crack surface films hinder hydrogen uptake to 
reduce da/dN and alter the proportions of each fatigue crack propagation mode. 


'Robert S. Piascik, formerly a graduate student in the Department of Materials Science at the University of Virginia, 
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I. INTRODUCTION 


A quantitative model for the kinetics of environmentally assisted fatigue crack propagation 
(FCP), based on mechanisms of crack tip damage, is central to alloy development and flaw 
tolerant life prediction 0 ' 2 *. Research reported in Part I characterized FCP in precipitation 
strengthened Al-Li-Cu alloy 2090 stressed in controlled inert (vacuum and helium), oxygen, and 
hydrogenous (pure water vapor, moist air and aqueous NaCl with varying polarization) environ- 
ments (3) . These experiments produced intrinsic fatigue crack growth rate (da/dN) versus stress 
intensity range (AK = K max - K min )data that are independent of crack closure and non steady- state 
chemical effects. A detailed fractographic analysis was presented in Part II; microscopic crack 
paths were defined and correlated with AK and the underlying microstructure (4) . Those 
environments that are capable of producing atomic hydrogen on clean crack surfaces enhanced 
rates of FCP in 2090, and produced unique intersubgranular and transgranular {100} 
crystallographic fatigue cracking modes (3,4) . These observations were qualitatively explained 
based on hydrogen environment embrittlement. 

The objective of the research in Part III is to understand and model the micromechanical 
and chemical bases for environmental FCP in Al-Li-Cu alloys. A framework is developed to 
model the relationship between da/dN and AK, based on crack tip damage related to the 
interaction between dissolved hydrogen, plastic deformation and microstructure. The model is 
used to explain the effects of important variables including AK, environment chemistry, 
microstructure and loading frequency. 

A. Background 

The hydrogen embrittlement mechanism to enhance fatigue crack propagation in aluminum 
alloys exposed to moist air, water vapor and chloride electrolytes is based on broadly compelling, 
albeit circumstantial, evidence^’ 5 ' 32 *. Atomic hydrogen is chemically adsorbed on clean crack 
surfaces as the straining aluminum alloy chemically reacts with gaseous H 2 0, or after electro- 
chemical cathodic reduction of hydrogen ions or water. Fatigue crack tip process zone damage 
is probably promoted by dissolved hydrogen through hydride formation, decohesion or localized 
plasticity based mechanisms (5,l7 ' 29,33,34) . Although evidence is lacking, environmental FCP 
mechanisms were proposed based on gaseous adsorption and crack surface film formation° 2,35 ' 37) . 
Adsorbed hydrogen or oxygen on the crack surface interact with strained atomic bonds to promote 
brittle fracture at lowered stress. Surface films may affect crack extension by inhibiting reversible 
slip, promoting slip homogenization, localizing the distribution and morphology of persistent slip 
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bands, or reducing crack tip welding. For alloys in aqueous electrolytes, a film forma- 
tion/rupture/dissolution sequence was modeled to predict da/dN-AK (12,38,39 \ however, this 
approach has not been applied to aluminum alloys. 

Quantitative models for hydrogen environment enhanced FCP kinetics must be based on 
one or more governing microscopic damage processes, and must integrate the mechanical and 
chemical components of crack tip damage. The mechanical contribution is described by relating 
the effective stress intensity to the continuum crack tip stress, plastic strain and strain rate fields' 

42) , and by determining the fracture process zone and damage sites. The chemical contribution 
requires modeling of occluded crack chemistry, crack surface reactions and hydrogen diffusion 
in the crack tip process zone (2) . Damage and da/dN depend on the amount of hydrogen delivered 
to the crack tip fracture process zone, and are described by a failure criterion between normal 
stress, strain and hydrogen concentration. Crack propagation occurs either per cycle or discon- 
tinuous^ after damage accumulation. To date, modeling has addressed either the crack chemistry 
of hydrogen production or crack tip plasticity damage without considering hydrogen; there is no 
integrated model for da/dN as a function of AK, chemical and microstructural variables. Results 
in each regard provide the basis for the current modeling. 

1 Mag§ Transport and Reaction Ra te Limited FCP 

Hydrogen embrittlement models were developed to predict da/dN, as simply proportional 
to the amount of atomic hydrogen delivered to the crack tip fracture process site per single load 
cycle. These models hypothesize that da/dN is rate limited by one of the following: (a) 
environmental mass transport in the crack environment, (b) crack surface chemical or electro- 
chemical reactions, or (c) hydrogen diffusion within the plastic zone. For the aluminum alloy- 
water vapor system, Wei and coworkers modeled da/dN based on slow mass transport within the 
occluded crack environment®. With limited exceptions, water vapor reactions on clean straining 
aluminum surfaces rapidly proceed to completion during a load cycle and are not rate limiting. 
Hydrogen production, and hence da/dN, are controlled by impeded molecular (Knudsen) flow 
within the crack, but only below a bulk environment "saturation pressure"' 28,43,44 *. Predictions 
of da/dN as a function of water vapor pressure, loading frequency and stress ratio (crack opening 
shape) are confirmed by experimental observations for several aluminum alloys, including 

2090 (3,28,43,45 ’ 48) 

The kinetics of hydrogen producing cathodic reactions on aluminum surfaces in aqueous 
electrolytes have not been characterized and related to environmental FCP (2,44) . For steels, Wei 
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and coworkers relate the corrosion fatigue component of da/dN to the amount of adsorbed atomic 
hydrogen produced per load cycle by cathodic H + and H 2 0 reduction reactions, coupled with 
transient anodic dissolution* 2 - 44 *. Hydrogen production and da/dN are equated to the ratio of the 
anodic charge transferred per loading cycle and the amount associated with completed reaction 
of the clean surface. Charge transfer is measured in a fracturing electrode experiment (49) . This 
approach reasonably predicts the loading frequency dependence of corrosion fatigue crack growth 
for ferritic steels in electrolytes* 2 ’ 44 - 49 *. The model has not been applied to the aluminum al- 
loy/aqueous chloride system where surface reaction rates are likely to be fast relative to typical 
loading frequencies* 45 *. 

Contrary to the surface reaction model, environmental FCP rates may be controlled by 
atomic hydrogen diffusion within the cyclic plastic zone* 6 - 50 *. Da/dN is assumed to be proportion- 
al to the per cycle penetration distance of hydrogen (Ax), and thus the square root of hydrogen 
diffusivity in aluminum and diffusion time (the reciprocal of the loading frequency). This 
relationship described the frequency dependence of the maximum crack growth rate where a 
hydrogen embrittlement induced intergranular to transgranular crack growth transition was 
observed* 6,50 *. 

Chemically based crack growth rate models can predict the dependencies of da/dN on 
environmental activity (electrode potential, gas pressure and temperature) and time (loading 
frequency). The effects of AK, mean stress and alloy microstructure are not described because 
da/dN is equated to hydrogen production, without considering process zone damage and the 
governing failure criterion. A critical assumption of the mass transport and reaction limited 
models is that the crack extends by a small increment per load cycle, however, there is no 
evidence to support continuous crack growth. Chemical models for hydrogen environment FCP 
do not address the uniquely damaging influence of cyclic deformation and the specific 
microscopic modes of process zone cracking. 

2. Process Zone Plasticity Damage During FCP 

The interaction of cyclic plastic strain with microstructure and dissolved hydrogen near 
the crack tip must be understood to model environmental FCP behavior. Hydrogen volume 
diffusion in aluminum alloys is slow at room temperature, on the order of 10 13 to 10 15 m 2 /sec* 19 
21,51 55) , and may not support per cycle FCP over distances in excess of several pm. Arguments 
to the contrary were not proven* 6 ’ 50 *. Atomic hydrogen, residing in the cores of mobile disl- 
ocations, may sweep ahead of the straining crack tip up to 10 4 times faster than that enabled by 
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volume diffusion 07 ' 26 ' 27 54,56 5X) - Enhanced grain boundary hydrogen diffusion was also proposed, 
although little evidence exists for this mechanism 051 ' 56 - 59 '. Boundary precipitates such as T, are 
potential hydrogen traps that block a continuous diffusion path for hydrogen. 

Process zone microstructure plays an important role in hydrogen embrittlement of 
aluminum alloys. During cyclic loading, dislocation cell structures and microvoids form as plastic 
strain accumulates and when microstructural barriers (grain and subgrain boundaries, precipitates, 
and dispersoids) interrupt the movement of dislocations 1 *- 60 641 . Hydrogen atoms may segregate 
("trap") at these fatigue induced sites 12226,56 - 58,65,66 '. Crack growth rate is enhanced when the 
crack path coincides with trap sites. Trap sites that shield the crack path from hydrogen 
segregation slow da/dN. Embrittlement of aluminum alloys is presumably enhanced by localized 
slip planarity that concentrates both hydrogen and stress 1 

Considering mechanical fatigue damage based on plastic deformation at a crack tip, the 
unresolved issues are the precise nature of damage and the continuity of crack advance, be it per 
cycle or after accumulation over a number of load cycles 168,69 '. For da/dN above 10 mm/cycle, 
mechanical FCP is viewed to continuously occur during a load cycle and by a mechanical plastic 
blunting process that leads to microscopic striations. Laird proposed a crack blunting and 
buckling process for crack advance and resharpening 171 ". Pelloux argues that continuous and 
striated fatigue crack growth in aluminum alloys is produced as alternating shear concentrates 
along {111} slip planes and symmetric about a resultant { 100} crack plane 171 '. Both ductile {111} 
and brittle {100} striations have been reported, however, the role of environment is unclear for 
each process 16 - 30,10 ' 28,71 - 77 '. 

For da/dN below 10 4 mm/cycle, Lankford and Davidson propose incremental crack 
extension by a two stage process involving crack tip strain accumulation, blunting and incremental 
crack extension over a number of loading cycles 18,68,78 '. Incremental FCP was observed in 7075 
and MA-87 aluminum alloys stressed at near threshold cyclic stress intensities (AK^) in either 
moist air or nitrogen. The embrittling effect of moisture was ascribed to hydrogen which 
increased the increment of crack growth, and reduced the necessary failure strain and therefore 
the number of accumulated cycles for crack extension. 

For discontinuous FCP, several models quantitatively predict da/dN-AK dependencies 
based on a microstructurally allowable length for slip 141,42,61,68,69,78 '. The number of cycles for 
each increment of crack advance is given by a Cotfin-Manson-type failure criterion based on 
crack tip plastic strain range. While consistent with da/dN-AK data for fatigue in embrittling 
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moist air, these models do not consider hydrogen effects and are not broadly tested with crack 
closure-free data for intrinsic FCP in an inert environment' 68 *. Growth rate models have not been 
developed based on damage by the interaction of hydrogen and cyclic plasticity. 

3, Hydrogen Embrittlement and Cyclic leading 

Time-cycle dependent environmental FCP occurs below the threshold stress intensity for 
monotonic load cracking in all aluminum alloy-water vapor and chloride systems including 
2090 (3,m5) . Accordingly, cyclic loading and cyclic plastic deformation are uniquely damaging, 
as qualitatively explained by several hypotheses' 5 *: 

1. Fatigue crack convective mixing enhances the transport of deleterious species from the 
bulk environment to the crack surface. Inhibiting species at the crack tip are similarly 
ejected. 

2. Loading and unloading produce a sufficiently high crack tip strain rate to enhance film 
rupture, hydrogen production and hydrogen entry efficiency. 

3. Cyclic plasticity continuously transports atomic hydrogen by the dislocation mechanism 
from the crack surface-environment interface to strong trapping sites within the 
embrittlement process zone. 

4- Cyclic deformation produces crack tip dislocation cell or slip band structures that are 
strong trap sites and are embrittled by segregated hydrogen. 

5. High normal stress is maintained throughout the crack tip field because of continued crack 
tip sharpening by cyclic unloading. 

No model has been developed to quantify these explanations and to predict da/dN as a 
function of AK, chemical and microstructural variables. The critical role of cyclic dislocation 
plasticity (items 3 and 4, above) is emphasized in the following model. 

HI , Process Zone Model For Hydrogen Environment Assisted FCP 

Since mass transport and reaction to produce hydrogen have been emphasized, it is 
important to consider the interaction of this hydrogen with the highly localized distributions of 
stress, strain, dislocation plasticity and microstructure in the crack tip process zone. 

Presumably, hydrogen enhanced FCPdiscontinuously occurs over a distance (Aa), defined 
as the embrittlement zone (d ez ), and after a number of elapsed load cycles (AN). The average 
crack growth rate is given by: 

da/dN = (Aa/AN) = d e2 /AN (1) 
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where: Aa = fatigue crack advance distance 

AN = load cycles required for crack advance; a function, g(Ae p , a N , [H] pz ), of Af p , cr N and [H] pz 
^cl “ embrittlement or d ama & e 201,(5 at *be crac ^ hp; a function, f(Ae p , cr N ), of Ae p and a N 

A tp = crack tip plastic strain range; a function ot r, the distance ahead of the crack tip 
a N = crack tip stress normal to embrittled microstructure; a function of r 
|H | p7 = process zone atomic hydrogen concentration; a function of crack surface hydrogen 
concentration and crack solution electrochemistry. 

This "damage accumulation" approach was employed to model mechanical FCP based on the 
coupling of the crack tip strain field (either continuum or experimentally defined) with a strain- 
based fatigue failure criterion such as the Coffin-Manson rule< 8 - 4, - 42 - 78 TO >. The physical bases and 
quantitative relationships for Aa and AN are, however, likely to differ for FCP in hydrogen 
producing environments. 

For FCP in Al-Li-Cu alloys, the role of crack tip cyclic plastic deformation is to enable 
dislocation transport of hydrogen to supply process zone embrittlement. For most growth rates, 
concentration gradient bulk diffusion is insufficient to sustain environmental FCP. For example, 
for a bulk hydrogen diffusivity (D H ) of 5 x 10 15 m 2 /sec, diffusion throughout the cyclic plastic 

zone requires times equivalent to 60 load cycles for AK = 2 MPa/m, 250 cycles for AK = 3 
MPa/m and 6,000 cycles for AK = 7 MPa/m. 2 For all but the lowest AK levels, these 

numbers of cycles exceed likely values of hydrogen reduced AN. Since a diffusivity of 10 10 to 
10 " m 2 /sec is typical of dislocation transport of hydrogen in aluminum, particularly when slip 
is planar and localized, hydrogen rapidly penetrates the cyclic plastic zone; if mobile dislocation 

activity is sufficient* 57 '. 

Atomic hydrogen is accumulated within the crack tip process zone by cyclic plastic strain, 
but brittle cracking is caused by the crack tip normal, or opening direction stress, <r N . For 
fatigue crack advance, a N must be sufficiently high over the plasticity (Ae p ) based distance where 
hydrogen has accumulated. This distance is defined as the embrittlement zone, d ez . This view 
is assumed because normal stresses typically control interfacial (eg. SGC) and cleavage (eg. 
{100} crystallographic) cracking modes. Normal stress controlled FCP competes with cyclic 


hhe calculated cycles for hydrogen diffusion are based on elapsed time, t. at a 
the diffusion distance, 2(D H t) l/2 , equated to the cyclic plastic zone diameter, d 
reasonably approximated by the Rice-lrwin result: 


loading frequency of 5 Hz and with 
cpz . As discussed in Part II, d cpz is 


d cp7 = ( 1 / 1 2ir)( AK/a ysc ) 2 


with the monotonic yield strength, a ys , employed in place of the cyclic value (40) . 
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plastic strain-based damage, depending on the amount of process zone hydrogen and the level of 
depends on applied AK, crack closure (if present), and material cyclic flow properties; 
while a N depends on K max and material tr ys and work hardening exponent, n. Analysis of 
Equation 1 to predict da/dN-AK is complicated because At p , <r N and [H] pz vary with distance 
ahead of the crack tip and with time during one or more load cycles, and the "f" and "g" 
functions are undefined. Material flow properties may depend on dissolved hydrogen. It is, 
none-the-less, informative to develop estimates for the crack advance distance and the number of 
cycles for discontinuous crack advance. 

A. Determination of Crack Advance Distance - Aa 

For hydrogen environment FCP, the crack advance increment (Aa = dj is determined 
by the distance ahead of the crack tip where cyclic plastic strain (Ae p ) is sufficient to generate 
dislocation activity and enhanced hydrogen concentration. For heterogeneous planar slip 
materials such as aluminum alloy 2090, the nature of the dislocation cell structure is unclear (60 ' 80) ; 
the dominant process is assumed to be hydrogen transport by mobile dislocations to precipitate 
interface and subgrain boundary trapping sites. 3 The process zone model is simplified by the 
assumption that dislocation activity is defined by the peak plastic strain within the process zone. 
Differences in damage accumulation, that is dislocation activity and hydrogen concentration, 
associated with the crack tip strain gradient are not considered* 41 * - 79 *. 

Since significant plasticity is confined within a fraction of the crack tip cyclic plastic zone, 
it is reasonable to speculate that crack advance occurs over a fraction a of d • 

’ ’ CpZ* 

Aa = d cz = ad cpz ( 2 ) 

Considering this plasticity condition for d cz , a is defined by the critical plastic strain range (Ae p 
crit ) above which mobile dislocation activity is sufficient to transport H to trap sites. (Alternately, 
a critical strain may be required to produce a sufficient dislocation substructure for hydrogen 
trapping .) a is determined by equating Ae p _ cri( to the continuum mechanics description of the 
(stationary, plane strain) crack tip plastic strain range, Ae p , for the condition where r equals d ez 
or equivalently, ad cpz . Since: 

A ep = W ( V /E)[(d, p/ /r) l/(n+l) - 1] (3)4 


3 

For other alloys such as steels, a stable dislocation substructure is progressively formed by cyclic deformation. 
Hydrogen trapping at this substructure and associated microvoids may be dominant. Dislocation transport is not 
required for hydrogen in steels because of the relatively high hydrogen diffusivity of I0 10 m 2 /sec (58,81) . 

4 

Note that this continuum-based plastic strain distribution predicts a singularity at the crack tip where r equals 0. 
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where: 


a . = cyclic yield strength 
yac 

E = elastic modulus 
r = distance ahead ot the crack tip 
n = work hardening exponent 
oj = constant 

solving for a yields: 

a = [{(Ae crit E)/a y J + 1] 1 ^ 

where: l/(n+l) = 1 for peak aged aluminum alloys 

j (40) 

The continuum mechanics result in Equation 3 was experimentally demonstrated to describe the 
plastic strain distribution to within 1 pm of the crack tip. a, least fo, Fe-3%Si<»“ The result 

has not been verified for planar slip A 1 -Li alloys. 

a depends on alloy microstructure, particularly if the cyclic deformation mode varies, and 

environment, if the brittle crack path changes. While this parameter cannot be unambiguously 
specified, various estimations discussed in Appendix A yield values of a between 0.4 and 0.7. 

Normal stresses are sufficiently high at all positions within the embrittlement zone and do 
not alter Aa determined by dislocation plasticity and Ae cri( . As shown in Appendix B, continuum 
analysis of the crack tip opening stress distribution within d cp/ indicates that <r N at maximum load 
is constant and high (equaling 3 to 4 times <r ys , depending on work hardening) over a distance of 
3 to 250 times d ez , for R values of 0. 1 and 0.9, respectively. As applied K is cyclically varied 
from K min to K^, the maximum value of <r N is unchanged, however, the extent of the stress field 
penetration within the process zone increases from a location very near to the crack tip (Equation 
B1 with K min 2 ) to the value given by Equations B1 and B2 in Appendix B. Unloading from 
to K min produces a compressive singularity, with near-crack tip compression depending on the 
extent of crack blunting and superimposed upon the tensile field. Dislocation motion may reverse 
in response to compressive yielding, however, this action does not alter the conclusion that tensile 
stresses within the plastic strain-based d„. are sufficient to enable hydrogen cracking; Aa is defined 
by £ p <ni through a. 

B. Determination of Failure Criterion - AN 

In the process zone framework, AN is the number of load cycles required for dislocation 

transport to achieve a critical level of process zone hydrogen, [H] pz<rit , for brittle crack advance 
over Aa Here, failure occurs as a result of high cr N . This view of AN is stress-based and is 
fundamentally different from strain-based fatigue damage^ 1 The challenge is to predict 
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AN as a function of <r N and [H] pz . [H] pz depends on Ae p , the absorbed crack surface hydrogen 
concentration ([H] ABS ) in equilibrium with the occluded crack environment, alloy trap micro- 
structure and fatigue dislocation trap structure. Hydrogen is delivered to process zone trap sites 
by dislocation motion during the loading portion of the fatigue cycle; [H] pz increases with load 
reversals until the critical failure criterion between <r N and [H] pz ^ ril is achieved. Incremental crack 
growth occurs and the embrittlement process is repeated. The maximum value of <r N is constant 
during load cycling. A finite number of load cycles is required to saturate relatively strong trap 
sites within the process zone because of limitations associated with "pickup" of atomic hydrogen 
at the crack tip surface-environment interface, with "deposition" of hydrogen as the mobile 
dislocation nears the trap, and because of the high hydrogen capacity of strong traps. 
Qualitatively, AN should decrease as the following factors increase; 

1 . Mobile dislocation density, which increases with increasing Ae p and hence increasing AK, 
and which depends on slip morphology; 

2. a N which is essentially constant with increasing AK, but which increases with increasing 
a y , and n; 

3. Hydrogen environment activity and uptake efficiency which increase [H] ABS ; 

4. Damaging trap sites (e.g., subgrain boundary T, precipitates, stable dislocation cells, 
localized slip bands and microvoids). An increased density of energetically beneficial 
traps increases AN; 

It is not possible to present quantitative relationships that define AN as a function of these factors. 
These relationships could be experimentally defined by low cycle fatigue-type experiments with 
notched and hydrogen precharged specimens. Such data have not been reported. 

The load cycles for discontinuous crack extension can be estimated from measured crack 
growth rates for alloy 7075 <3) and values of d cz estimated through Equation 1. Figure 1 is a plot 
of this process zone model estimate of AN compared with results calculated by a measurement 
and modeling approach applied to FCP in alloy 7075 in inert and hydrogen producing environ- 
ments (8) . In this latter method, AN is calculated from in situ SEM measurements of crack tip 
strain on the specimen free surface coupled with a Coffin-Manson criterion for the effect of plastic 
strain range on cycles to failure. The constants in the failure criterion were determined from 
da/dN-AK data. AN was not directly measured by the SEM method. For the process zone 
estimates, a values of 0.7 and 1 were used for moist air and vacuum, respectively. 

The process zone model estimates for AN shown in Figure 1 are in good agreement with 
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the SEM-based results 18,10 '. AN exceeds one, indicating discontinuous crack advance and AN 
generally decreases with increasing AK. The model and SEM results show fewer load cycles per 
crack advance event for moist air, an embrittling environment, compared to the inert environ- 
ment. (Da/dN for nitrogen containing 10 ppm H 2 0 are intermediate because low contaminant 
levels of H 2 0 are embrittling, as indicated by accelerated fatigue crack growth rates* 3,9,4 ,8 '.) 
Consistent with the previous discussion of slow hydrogen transport in aluminum, the small AN 
values in Figure 1 suggest that bulk process zone hydrogen diffusion must be augmented by 
mobile dislocations to sustain rapid environmental FCP rates. 

The specific AN versus AK relationship depends on the process zone failure criterion, for 
hydrogen producing moist air, damage is governed by [H] pz accumulation and a N -based cracking, 
while FCP in vacuum progresses by a plastic strain-based damage process* 42 '. Specific functional 
relationships between AN and Ae p cnl are proposed in an ensuing section which discusses the 
relationship between da/dN and AK. 


IV. Discussion 

A. Summary of Environmental FCP in Al-Li-Cu Alloy 2090 

Table 1 summarizes the microscopic crack paths, and the associated da/dN-AK power law 
relationships for alloy 2090 in various environments, as reported in Parts I and H (3,4,45) . The 
{111}, subgrain boundary (SGC) and crystallographic {100} cracking modes are listed by 
increasing environmental susceptibility; vacuum, helium and oxygen, all produce equally slow 
da/dN; followed by moist air, water vapor and anodic NaCl which promote the faster FCP rates. 
For nonhydrogenous vacuum, helium and oxygen, a single {111} cracking mechanism is 
observed and da/dN is characterized by a single AK 4 0 power law. For hydrogen producing moist 
air, water vapor and aqueous NaCl, changes in the da/dN-AK relationship are associated with 
different proportions of brittle SGC and {100} cracking. At low AK, where a AK 5 1 power law 
is observed for both anodic NaCl and moist air, only transgranular {100} cracking is observed. 
SGC occurs at moderate to high AK; this fracture mode is described by a single slope power law 
of AK 2 7 for anodic NaCl and AK 2 4 for water vapor. Slip plane cracking is produced at high AK 
for the aggressive environments. Complex changes in the da/dN-AK slope correlate with multiple 
fracture modes. For example, the da/dN-AK transitions and apparent "plateau" observed for 2090 
in moist air at moderate to high AK are associated with varying proportions of {100}, SGC and 
{111} slip plane cracking. 
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Table 1: Environmental Fatigue Cracking in Alloy 2090**’ 4 * 


Environment 

Region A 
Low AK 
( £ 2 MPav^ m) 

Region B 
Moderate AK 
(3 to 7 M Pic/ in) 

Region C 
High AK 

(> IOMPa\/m) Comments 

Vacuum 

{no 

{in} 

{no 

AK 40 

Helium 

{in} 

{in} 

{in} 

AK 40 

Oxygen 

{in} 

{in} 

{in} 

AK 40 

Moist Air 

{100} 

{ I00}/SGC/{ 111} 

SGC/{ 111} 2 

AK' 1 (low AK) 
AK 1 5 (mod AK) 
AK 3 ' 7 (high AK) 

Water Vapor 

{100} 

SGC 

SGC/{111} 3 

AK 2 4 (mod AK) 4 
AK 2 4 (high AK) 4 

NaCI (anodic) 

{100} 

SGC 

SGC/{111} 3 

AK 5 1 (low AK) 
AK 2 7 (mod AK) 4 
AK 2 7 (high AK) 4 


1 Defined in Figure 13, Part II <4 \ 

2 Increasing proportions of mechanical {111} cracking are observed with increasing AK. 

3 Primarily SGC with small proportions of {111} cracking at high AK. 

4 Single slope behavior based on a single fracture mode. 

B. The Effects of Mechanical, Microstructural and Chemical Variables 

Considering the environmental FCP behavior of alloy 2090 described in Parts I and II and 
summarized in Table I (3,4,45) , the hydrogen embrittlement process zone framework explains the 
measured AK and environment activity effects on the microscopic crack path, and the complex 
relationship between da/dN and AK. 

1. Why Poes the Intrinsic Fatigue Cracking Mode Depend on AK? 

A key aspect of environmental FCP in Al-Li-Cu and other 2000 series aluminum alloys 
is the complex dependence of crack growth rate on AK, but only for aggressive environ- 
ments <3,4,45) . Mechanical explanations for this behavior; including crack closure, static mode 
cracking above K ISCC , stress state, and slip barriers; are not important* 3 *. Rather, the relationship 
between da/dN and AK is governed by AK-induced transitions in the microscopic modes of 
hydrogen enhanced FCP. 

The results in Table 1 (and in Figures 12 and 13 of Part II) establish that the microscopic 
mode transitions between Regions A, B and C are controlled by the interaction between the 
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embrittlement zone volume and the included microstructure. Unlike previous arguments based 
on slip damage of unspecified form over a critical distance* ’ ’ \ the process zone model 
suggests that microstructure controls hydrogen segregation by trapping and associated 
microcracking. Figure 2a is a scale drawing of the relationship between the L-T crack tip process 
zone (d ez = ad cp/ ) and the 2090 subgrain structure at high AK (Region C). The average subgrain 
size of 5 /im (transverse direction) is shown. At a AK of 15 MPav/ m (K max = 17 MPa\/m), the 
embrittlement zone volume (d cz = 17 /tm) encompasses many subgrains in the crack growth 
direction and along the crack front. During cyclic loading, atomic hydrogen produced at the 
crack tip surface is transported by mobile dislocations to strong trap sites at subgrain boundaries 
within this large embrittlement zone. Bulk hydrogen diffusion is minimal at high AK because few 
load cycles and a short time are required to hydrogenate the process zone by dislocation transport 
and to achieve the local failure criterion for discontinuous crack extension. In this regime 
mechanical slip band cracking competes with environmental SGC. Aggressive environments that 
produce high levels of crack surface hydrogen favor brittle intersubgranular cracking and little 
or no SBC. For less aggressive environments such as moist air, some SBC occurs, perhaps 
stimulated by hydrogen* 3 '• 57) and in parallel with SGC. 

For near-threshold environmental FCP (Region A), embrittlement of subgrain boundaries 
is not possible because such boundaries are not hydrogenated by the localized dislocation activity 
associated with the small cyclic plastic strain-based process zone; Figure 2b. At these AK levels, 
d ez is not sufficiently large for dislocation transport of hydrogen to reach subgrain boundaries. 
In the absence of SGC, environmental FCP proceeds along {100} planes. As AK increases, the 
embrittlement zone encompasses the smallest subgrains. The boundaries are enriched with 
hydrogen by trapping and the proportion of SGC increases while {100} cracking within subgrains 
decreases. The growth rate transition (A to B) and the beginning of SGC are expected at the AK 
level where d cz just equals the smallest subgrain size. Within Region B, {100} and SGC occur 
in parallel in proportions that are governed by the relationship between the distribution of 
subgrain sizes and the extent of dislocation transport within d cz . 

For Region A, a complication exists because both dislocation transport and bulk hydrogen 
diffusion may contribute to crack tip embrittlement due to the large number of load cycles, and 
the associated time, for each crack advance event. For example, the process zone model suggests 
that on the order of 200 load cycles are required to achieve the hydrogen-normal stress failure 
criterion and incremental crack growth over a relatively small advance distance (see Figure 1). 


13 



For a diffusivity of 10 14 m 2 /sec, the diffusional penetration of hydrogen per increment of crack 
growth (e.g., 200 load cycles at a AK of 2 MPa\/m and 5 Hz) is 1 n m. D H is uncertain and 
hydrogen diffusion may be slowed within d iv by dislocation and microvoid traps in the highly 
strained plastic zone <58) . The {100} to SGC mode transition begins at that AK level where 
dislocation transport occurs over a distance of the smallest subgrain diameter; but only for those 
cases where hydrogen supply by bulk diffusion is nil. This case is likely for high loading 
frequencies or fast crack growth rates (low elapsed time per crack advance event), and for low 
concentrations of crack surface hydrogen. When bulk diffusion of hydrogen is substantial within 
the process zone, intersubgranular cracking persists to AK levels below the transition value from 
the plasticity consideration. The trends in Figure 13 (Part II) are consistent with these arguments. 
Fora weak hydrogen environment such as moist air, SGC begins when d„ equals the smallest 
subgrain size. SGC is observed for AK and process zone size conditions that are smaller than 
the subgrain size for the aggressive NaCl environment. Above this transition, but for process 
zone sizes less than the subgrain size, dislocation transport augments bulk hydrogen diffusion, 
perhaps by translating the surface hydrogen source within the process zone and closer to, but not 
reaching, subgrain boundaries. The transition to { 100} cracking and the Region A da/dN-AK 
dependence indicate that bulk diffusion is insufficient to supply subgrain boundaries below a 
transition AK level and process zone size. Since monotonic load environmental cracking does 
not occur at the stress intensities where FCP is affected by the bulk hydrogen diffusion 
mechanism, cyclic loading is critical. In the absence of dislocation transport control, cyclic 
deformation may be required to destabilize crack surface protective films, thus promoting cathodic 
hydrogen production and entry efficiency. 

2. Why Does Hydrogen Promote flOOl and Intersubgranular Cracking? 

It is only possible to speculate on the origin of hydrogen enhanced {100} and intersubgran- 
ular cracking in alloy 2090; both microstructural features provide for local hydrogen trapping and 
therefore preferential fracture. Hydrogen trapping in complex aluminum alloys is poorly 
understood compared to steels (58 86 88) . High angle grain boundaries; dislocation structures and 
microvoids from fatigue damage; incoherent interfaces between the aluminum matrix and 
constituents, dispersoids or overaged precipitates; and somewhat more coherent interfaces 
associated with segments of strengthening and subgrain boundary precipitates are likely sites for 
irreversible hydrogen trapping in Al-Li-Cu-Zr alloys. Coherent interfaces, typical of unrecrystall- 
ized subgrains, under to peak aged intragranular precipitates, and isolated dislocation cores are 
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weaker reversible trap sites*™' 86 '. These later sites contain hydrogen in equilibrium with lattice 

dissolved solute and will yield hydrogen to stronger traps. 

Considering the SGC mode, it is reasonable to speculate that low misorientation angle 
subgrain boundaries decorated withT, precipitates* 4 ' provide effective sites for hydrogen trapping. 
The strength of the subboundary-T, trap cannot be assessed, however, it is likely to be of 
moderate strength* 32 ' 87 '. Hydrogen may be trapped at the incoherent ends of T, plates along 
subboundaries. The subgrain boundary may trap hydrogen, without T,, because of interfacial 
disorder or lithium solute segregation. As fatigue plasticity and damage accumulate at subgrain 
boundaries and T„ including within the soft 8 ' precipitate free zone, the strength of trap sites 
increases and perhaps approaches irreversibility. Subgrain boundary T, is a stronger trap site 
than intragranular T, and 8 ', because of the more disordered interface and localized fatigue 
damage, and is stronger than dislocation trap sites* 8 *'. Mobile dislocations will therefore transport 
hydrogen from the crack-environment interface, through precipitation hardened aluminum grains, 
to subgrain boundaries where the solute is deposited at the stronger trap sites. FCP occurs along 
the subgrain boundary path because the localized failure criterion between hydrogen concentration 
and normal stress is achieved at such sites because of hydrogen segregation. 

Subgrain boundary T, is likely to anodically dissolve in an aluminum matrix, and to 
concomitantly stimulate local atomic hydrogen, entry when such boundaries intersect the crack 
surface and encounter electrolyte* 32 ' 89 '. This behavior could promote hydrogen enhanced SGC, 
however, it is not the governing factor because similar intersubgranular cracking is observed for 
2090 in aqueous NaCl, very low pressure water vapor and moist air* 3 ’ 4 '. The absence of an 
electrolyte for the gaseous environments precludes T, dissolution, but not hydrogen trapping. 

The reason for hydrogen enhanced cracking along {100} planes is more problematic. This 
crystallographic fracture mode has been reported for stress corrosion cracking and corrosion 
fatigue in a variety of 2000 and 7000 series aluminum alloys (Section I, Part II)* 4 ' and for 
constant load liquid metal embrittlement of Al-Cu alloys* 90 '. Meletis recently reported 
transgranular "cleavage" cracking for hydrogen precharged Al-Li-Cu, without specifying the 
crystallographic crack plane* 32 '. Brittle cracking along {100} habit planes, typical of Guinier 
Preston zones or 0’ precipitates in Al-Cu alloys, does not explain {100} cracking in 2000 and 
7000 series alloys that contain precipitates with either {111} habit planes (e.g., T, or fi) or that 
are spherical without a well defined planar interface with the aluminum matrix (e.g., 5’)* 90 '. A 
lattice decohesion mechanism provides a rational, but speculative, basis for {100} cracking. 
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Hydrogen may accumulate in octahedral lattice sites that lie within { 100} planes in the aluminum 
lattice. The source of this hydrogen is short range diffusion (or "leakage") from the trapped 
accumulation on {111} slip planes in planar slip alloys such as 2090. Crystallographic cracking 
results from crack tip normal stresses acting on disrupted lattice bonding about {100} planes. 
Atomic bonding energies in aluminum are reduced by solute hydrogen (l7) . Both hydrogen 
accumulation and local tensile stresses could be enhanced by a dislocation mechanism based on 
interacting {111} glide planes to produce unresolvable microvoid damage nominally oriented 
along {100} planes* 7 1 ,9 l9J) . 

The recent report of the LiAlH 4 phase, formed in alloy 2090 during cathodic polarization 
(-1500 mV M in 3% NaCl with a hydrogen recombination poison) and kinetically stable at 23°C, 
provides a speculative explanation for environmental FCP by the {100} and SGC modes <94,95) . 
Submicron sized hydride plates appear to preferentially nucleate at 5 (AlLi) precipitates, and by 
inference at T, precipitates, both within subgrains and at subgrain boundaries. The hydride is 
thermodynamically stable compared to common lithium bearing precipitates in Al-Li-Cu alloys. 
Subgrain boundary hydride provides a logical reason for preferred hydrogen localization and SGC 
within the crack tip process zone. Intragranular LiAlH 4 forms as plates on {100} habit planes 
in the aluminum matrix. TEM observations reveal that hydrides crack along {120} planes in the 
hydride and nearly parallel to {100} planes in Al (94 95) . No fractographic evidence was obtained 
during the current study to support a hydride-based mechanism and the original study did not 
demonstrate hydrides on the intersubgranular or transgranular "cleavage" surfaces produced by 
stress corrosion cracking in the 2090/NaCl system. None-the-less, extensive evidence 
demonstrates that hydrides are likely to form in AI-Li-Cu alloys <94 ’ 95) . The role of this phase in 
environmental FCP must be further evaluated. 

The transgranular {100} cracking mode transitions to intersubgranular cracking because 
of the strength of subboundary-T, hydrogen traps. As AK increases so that the process zone just 
encompasses a substantial number of subgrains, hydrogen associated with the {111} slip plane 
and {100} trap sites redistributes to the stronger sites at subboundaries. Without hydrogen 
accumulation, {100} (and {111}, if operative) cracking does not occur. Without hydrogen 
redistribution to strong subboundary traps, {100} cracking would persist to high AK levels, in 
spite of the fact that the process zone intersects subgrain and grain boundaries. The transition 
from crystallographic to interface hydrogen cracking in aluminum alloys depends on both AK (d ez ) 
and hydrogen partitioning between trap sites. Temperature will affect hydrogen trapping and the 
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operative crack path. 

3. Whv Does the Hydrogen Environment Affect the Functional Relationship Betw een 
ria/dN and AK? 

The relationship between da/dN and AK for a hydrogen-producing environment is 
determined by the AK dependencies of Aa, AN and the proportions of the { 100}, intersubgranular 
and {111} cracking modes. For each environment, one or more power-law relationships 
superimpose to define the crack growth kinetics law. 

Since Aa equals ad^, based on dislocation activity control of crack advance (Equation 1), 
Aa should increase with AK 2 , independent of the microscopic mode of hydrogen cracking. 
Specifically, from Equations 2 and 4, and d cp/ in Footnote 2: 

Aa = AK 2 [{Af cril E + a >s }<" + 'V^xa^" + 3) ] (5) 

Aa is determined only by the extent of mechanical slip within the process zone and is independent 

of hydrogen concentration and tensile stress, u N . 

AN decreases with increasing crack tip dislocation activity that provides increased 
hydrogen uptake to the embrittlement zone, per load cycle and for a constant crack surface 
environmental condition. That is, the number of load reversals required to achieve the critical 
trapped hydrogen concentration for crack advance decreases as the dislocation density within the 
embrittlement zone increases. Based on this argument, it is reasonable to assume that AN 
inversely depends on the total number of mobile dislocations within the crack tip embrittlement 
zone volume, N m , raised to some power, -<t>. Assuming that mobile dislocation density is propor- 
tional to plastic strain <96) , N m is described by the integral of the plastic strain range distribution 
from the crack tip to d c/ . Phis calculation, summarized in Appendix C, shows that: 

AN oc AK 2 * (6) 

This approach to estimating the AK dependence of AN is supported by the observation that 
intrinsic rates of environmental FCP generally depend on AK raised to a single power, at least 
for specific regimes of crack growth. The maximum crack opening tensile stress (<r N ) in a load 
cycle depends on <r ys and N, but not on AK; the stress aspect of the hydrogen failure criterion 
does not introduce an additional AK-dependent term. AN depends on the activity of the hydrogen 
environment. Combining the AK dependencies of Aa and AN yields: 

Aa/AN oc AK 2 /AK ^ oc AK <2 + 2<w (7) 

Equation 6 describes the AK dependence of da/dN for a single microscopic damage 
process, however, macroscopic FCP rates can be determined by multiple concurrent damage 
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mechanisms. Wei and coworkers argue that the total environmental fatigue crack growth rate 
(da/dN ENV ) equals the sum of the rates for individual mechanically and mechanical-chemically 
driven processes, that operate in parallel and each uniquely dependent on AK, times the areal 
fractions (9j) of each process (97) . For the case of alloy 2090 in the various environments: 

da/dN ENV oc e {l | l) AK 2/tf + 0 {m AK« + 2 *< ,00 » + 0 SCC AK (2 + 2 ^ SGC > (8) 

The first term describes mechanical FCP as typified by cyclic loading in vacuum. For this case, 
damage accumulation modeling predicts a power-law exponent of2//3, where 2 reflects the AK- 
dependence of the crack tip plastic strain range (Equation 3) for a nonhardening alloy, is the 
Coffin-Manson exponent employed as the failure criterion relating AN to Ae p , and crack advance 
occurs over a constant distance independent of aK < 4I HJ,84) . The AK-dependent process zone 
perspective applied to purely mechanical fatigue predicts that AK is raised to the power (2 + 
2//3). /3 has not been predicted from a damage model, but is measured to equal 0.5 +. 0.1 for 
alloys. Accordingly, FCP in a truly inert environment should depend on AK raised to between 
the fourth and sixth powers, depending on the AK-dependence of the damage distance; an 
exponent of 4.0 is observed for alloy 2090 (Table 1) (3) . The second and third terms in Equation 
8 describe {100} and intersubgranular cracking, respectively, with each uniquely depending on 
AK through Equation 7. 

When a single microscopic mode dominates FCP fora range of AK values, as is the case 
for 2090 in NaCl with anodic polarization (Figures 12 and 13, Part II), the amounts (9;) of the 
minority modes approach zero and da/dN depends on a simple power-law function of AK. This 
case is idealized in Figures 3a and 3b for two constant levels of crack surface absorbed hydrogen 
concentration, [H] ABS1 and [H] ABS2 . Since Regions A and B/C correspond to predominantly {100} 
and SGC cracking, respectively, the 4> component of the power law exponent in Equations 7 and 
8 must depend on the microscopic mode of cracking, ^-values for {100} and SGC cracking (1.6 
and 0.4, respectively) are calculated from the measured slopes of Region A and B/C (5. 1 and 2.7, 
respectively, Table 1). 

<t> may depend on the microscopic cracking mode for two reasons. A constant amount of 
dislocation activity (N m ) could differently supply hydrogen to the subgrain boundary and {100} 
fracture sites. That is, the efficiency of hydrogen transfer from dislocations to hydrogen traps 
at lattice and subboundary locations may vary. Hydrogen transfer between trap sites is likely to 
depend on the relative hydrogen binding energies for each site (58,87) , and on the dislocation 
structure on the active slip plane, which should vary with dislocation density and applied AK. 
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Secondly, the critical concentration of hydrogen required for brittle crack advance, at a given <r N , 
will vary for the two fracture mechanisms. The relationship between AN and N m will vary to 
reflect these physical processes. 

4, Why Does \H1 ab s Affect t h t> Magnitude and AK Dependence of da/dN ? 
Comparison of FCP in NaCl (anodic polarization) to cracking in the moist gas and cathodic 
polarization chloride environments shows that decreasing environmental production of crack 
surface hydrogen results in reduced absolute crack growth rates, a more complex da/dN-AK 
dependence with "plateau" behavior, and increased transition AK levels. This behavior is 
explained within the process zone framework and the assumption (Section B. 6) that decreasing 
concentrations of absorbed atomic hydrogen are produced by A1 surface reactions with NaCl 
(anodic polarization), water vapor, moist air and NaCl (cathodic polarization). 

The number of load cycles required to achieve the critical level of process zone hydrogen, 
at a given AK and hence mobile dislocation density, decreases with increasing crack surface 
absorbed hydrogen concentration in equilibrium with the local environment ([H ] ABS ). This 
hypothesis is based on the fact that the hydrogen capacity of strong traps within the process zone 
is independent of [H] ABS ; with increasing [H1 ABS . decreased surface-to-trap site activity of 
dislocations (NJ is required to produce [H] cn , Crudely, the product of [H] ABS and N m is 
constant. The increment of fatigue crack advance, Aa, is not affected by varying [H] ABS . Figures 
3a and 3b illustrate the simple effect of decreasing [H] ABS on da/dN-AK when a single 
microscopic fracture mode operates. For FCP governed by dislocation transport of hydrogen 
(Figure 3a), decreasing [H] ABS (from [H] ABS1 to [H] ABS2 ) depresses each power-law segment to 

lower da/dN at any AK level; the transition AK between {100} (0, = 0 {1O o} = L0 > and SGC (02 
= f) = 1.0) is not changed. Bulk diffusion may augment hydrogen transport during FCP at 

near threshold rates, as discussed in Section B. 1 . For the higher hydrogen activity environment, 
the dislocation and bulk hydrogen transport-based transition AK shifts to lower values and the 
power-law segments increase at constant slope with increasing hydrogen uptake; Figure 3b. 

C n. Multifile Fmctnrp Modes Affect The da/dN-AK Relationshlg? 

For some hydrogen environments, two or more damage modes can operate in parallel over 
a range of AK; Equation 8 must be evaluated in tom and complex da/dN-AK relationships are 
likely. Fractographic results establish that multiple cracking modes occur when the hydrogen 
environment activity is reduced compared to anodic aqueous chloride. This result is reasonable 
because as AN increases with decreasing [H] ABS at any AK, both SGC based on bulk hydrogen 
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diffusional supply and SBC based on mechanical damage accumulation assume increased 
importance. For multiple-mode FCP, the da/dN-AK dependence suggested by Equation 8 is 
represented by the bold lines in Figure 3c. liach of the three possible microscopic cracking 
modes ({ 100}, SGC and {111}) for alloy 2090 are associated with a single power-law da/dN-AK 
relationship, as indicated by the dashed lines for unity values of 0,, 0 2 and 0 3 , respectively. As 
AK increases, {100} cracking coexists with an increasing proportion of SGC; the two power-laws 
are added including increasing 0 2 and decreasing O,. Since {100} cracking is characterized by 
a higher power-law exponent, the sum of the rates for the parallel processes may exhibit a 
maximum. When only SGC occurs (0 2 ), da/dN-AK follows the single power law indicated by 
the short segment in Figure 3c. With further increases in AK, {111} SBC occurs and replaces 
SGC in increasing proportions. Since slip band cracking occurs at slower rates relative to 
intersubgranular FCP, Equation 8 indicates that da/dN ENV follows one of the dashed bold lines, 
depending on the relative values of 0 2 and O,. A single power-law da/dN-AK relationship is 
produced when 0 3 approaches unity, or when rates of SGC are inconsequential compared to slip 
band fatigue cracking. The concepts embodied in Equation 8 demonstrate the basis for complex 
da/dN-AK relationships, including "plateau" behavior. 

An example prediction of multiple fracture mode da/dN-AK "plateau" behavior is 
demonstrated in Figure 4 for peak aged alloy 2090 in moist air. Here, measured (1) and predicted 
(Equation 8) intrinsic FCP rates are in excellent agreement. Predicted crack growth rate 
(da/dN ENV ) is written in terms of the observed fatigue crack growth modes according to Equation 
8 : 

da/dN Moisl Air OC 0{| 1 1 j (da/d N| (Jorl ) -I- 0{ |(X)}(da/dN{ 1(X)}N„('|) 3" ®SGc(^ a /dN Water Vapor) (9) 

The second and third terms are empirically known because cracking in the 2090 was purely {100} 
for NaCl at low AK and purely intersubgranular for water vapor at intermediate AK. 5 
Substituting the measured da/dN-AK power-law exponent (Table 1) for each fatigue crack growth 
mode gives: 

da/dN MoirtAir = 0 {m} (3xlO 9 )AK 40 + 0 {1(X)} (lxlO*)AK 5 ' + 0 SG c(2xlO- 7 )AK 2 4 (10) 

with AK in units of MPa\/ m and da/dN in mm/cycle. In contrast to the exponents, the power-law 
coefficients decrease with decreasing [H] ABS (Figure 3a) and are not independently known. The 

values in Equation 10 were determined by fitting the portions of the moist air prediction to the 


5 lt is equally valid to describe SGC by intermediate AK FCP observed for NaCl, Table 1. 
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data in the AK regimes where single cracking modes dominated. 0, values were computed by 
a trial and error comparison of the data and model prediction in Figure 4. Figure 5 contains the 
distribution of 0 values for {111}, {100} and SGC modes used to calculate da/dN Moist Aif in 
Figure 4. These distributions are reasonable based on tractographic estimates of 0 s for Regions 

A, B and C. 

Quantitative fractography is required to verify the predicted fractions of each cracking 
mode shown in Figure 5. particularly in Region B where the proportions of each mode are 
complex and not easily discernable. The prediction in Figure 5 shows that complex da/dN-AK 
relationships are caused by varying proportions of microscopic cracking modes. While the power 
law exponent for each mode appears to be independent of hydrogen environment activity, the 
preexponential coefficients are not. Additional research is required to independently predict the 
effect of (H] abs on AN, and hence the absolute values of da/dN at any AK. 

Multi-slope da/dN-AK and plateau crack growth rates are reported for environmental FCP 
in a variety of ferrous and 7000 senes aluminum alloy-aqueous chloride systems 01 . Commonly, 
the transition AK is associated with a change in the fracture mechanism, and the governing factor 
is argued to be mechanical (or stress intensity) rate limitation at low AK and chemical (or 
hydrogen production and diffusion) limitation causing the plateau at higher aK , 6O0 '"' ,,) . This 
notion is supported by the frequency independence of low AK growth rates and by the strong 
frequency dependence of da/dN within the plateau for 7000 series aluminum alloys and steels 0 ’. 
For 2000-series alloys 0 ’, including alloy 2090°’, in aqueous chloride and saturation pressure 
water vapor, da/dN does not increase with decreasing loading frequency (see Section IV.B.b.b). 
Accordingly, the complex da/dN-AK relationships depicted in here are not explained by chemical 
rate limitations. While the shape of the da/dN-AK law was ascribed to transitions in fracture 
mode, detailed measurements and models were not proposed' . The process zone framework, 
Equation 8, provides a rational explanation. When operative, rate limited environmental FCP 
further complicates the intrinsic da/dN-AK relationships predicted by the process zone model for 
multiple parallel cracking modes®. Additionally, extrinsic crack closure will complicate da/dN- 

AK laws for both inert and embrittling environments . 

6 rv~ .h. tw Zone F r ~™-wnrk Explain the Effects of Hydrogen Environ ment 

Artivitv and Loading Fre quency? 

Data in Parts 1 and II and the process zone damage framework suggest that decreasing 
concentrations of absorbed hydrogen on alloy 2090 crack surfaces are produced in order by NaCl 
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(anodic polarization), pure water vapor, moist air and NaCI (cathodic polarization). Da/dN, the 
proportions of each microscopic cracking mode, and the da/dN-AK dependence depend on 
[HJabs- For th e moist gaseous environments, molecular transport modeling and experimentation 
establish that [H] ABS and hence da/dN increase with: a) increasing exposure (water vapor pres- 
sure/loading frequency), but only up to a saturation level, b) increasing R, and c) decreasing 
oxygen partial pressure and blocking surface oxide' 2 - 43 - 46 - 48 ’. These expectations are supported 
by FCP data for alloy 2090' 3 - 45 ’ and are consistent with the process zone concepts of dislocation 
transport and microstructural trapping of hydrogen. Process zone hydrogen diffusion does not 
limit FCP because of dislocation enhanced mass transport and the short diffusion distance over 

A ( 6 , 50 ) 
u ez 

Crack growth in aqueous chloride is qualitatively understood based on occluded crack 
electrochemistry and crack surface passive films that reduce hydrogen entry to the process zone. 

a. Electrode Potential: For alloy 2090 in aqueous NaCI, the important result to explain 
is the substantial rates of FCP under anodic polarization compared to the inhibiting and crack 
arresting effects of cathodic electrode potentials' 3 ’. 

A schematic of crack tip electrochemical reactions is shown in Figure 6a. For anodic 
potentials near the free corrosion value, dissolution occurs under relatively surface film-free 
conditions and the resulting cations (predominantly aluminum) are hydrolyzed to decrease crack 
tip pH and increase hydrogen production by proton reduction' 100 - 102 ’. As potential decreases, 
dissolution rates decrease, Al +3 hydrolysis is reduced, pH increases and [H] ABS decreases, as 
indicated in Figure 6b. With decreasing potential, water reduction becomes the dominant source 
of hydrogen. The minimum in [H] ABS correlates with monotonic load hydrogen cracking in 
aluminum alloys at both anodic and cathodic potentials' 51,95 ’. 

The complicating effect of a crack surface film must be considered for Al-Li alloys. At 
cathodic potentials, the increasingly alkaline crack solution produces a black surface film, possibly 
Al(OH) 3 and LiOH' 3 ’. Speculatively, LiOH hinders atomic hydrogen permeation into the process 
zone; [H ABS ] is lowered, thus explaining the beneficial effect of cathodic polarization. While 
reduced da/dN is explained based on either the chemistry argument represented in Figure 6, or 
by the film effect, the prolonged times for fatigue crack growth reinitiation when changing from 
cathodic to anodic potential' 3 ’ are evidence for the dominant effect of a protective crack tip surface 
film. 

Within the context of the process zone model, the important conclusion is that cathodic 
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polarization reduces [H] ABS for alloy 2090 in NaCI. Decreasing crack surface hydrogen 
concentration results in "plateau" behavior because: a) the magnitude of da/dN is reduced for both 
pure { 100} and intersubgranular cracking, b) the {100} to SBC transition AK increases because 
of decreasing bulk diffusion of hydrogen, and c) SBC is increasingly important at higher AK 
when rates of SGC are reduced. This behavior is illustrated by the bold curve in Figure 3c. 

b. Loading Frequency: Considering alloy 2090 in NaCI (anodic polarization), the 

important result to explain is that loading frequency (f) has no effect on da/dN for near threshold 
AK levels and that da/dN moderately increases with increasing frequency for higher AK ( '. 

The frequency dependence of da/dN for environmental FCP in aluminum alloys is 
determined by the kinetics of crack environment mass transport, chemical or electrochemical 
surface reactions and process zone hydrogen diffusion, and by crack tip strain rate effects on 
surface reactions, film stability and [H1 ABS ' 2,) . The measured frequency dependencies of da/dN 
for alloy 2090 in aqueous chloride suggest that: (a) surface anodic dissolution and cathodic H 
production reactions progress to completion for the frequences considered (0.07 Hz to 20 Hz) 
and (b) process zone hydrogen diffusion is fast during the multiple load cycles necessary for crack 
advance. Accordingly, environmental FCP in the 2090/NaCI system is not transport or surface 
reaction rate limited. Rather, we speculate that increased crack tip strain rate causes increased 
hydrogen production and/or uptake to reduce AN and increase da/dN. 

Fast surface reaction kinetics are consistent with results reported by Wei and coworkers 
for 7000 series aluminum alloys' 2 - 28 ' 44 '. The process zone model suggests that hydrogen transport 
ahead of the crack tip is not rate limiting because of rapid dislocation transport and short diffusion 
distances, d cz . For hydrogen transport in the process zone to be growth rate limiting, a^ loading 
frequency of 150 Hz or higher is required; based on a hydrogen diffusivity of 10 m /sec for 

dislocation transport, a process zone size of 7 ,m (AK = 10 MPa/m) and a AN of 20 cycles. 
At low AK, diffusion distances are smaller, AN is larger, and extremely high loading frequencies 

(> 150 Hz) are required to limit process zone hydrogen diffusion. Stanzl etal. demonstrated that 
relatively high pressure pure water vapor accelerates low AK FCP rates in 2024-T3, compared 
to vacuum, and at an ultrasonic loading frequency of 20,000 Hz' 103 '. A complex da/dN-AK 

relationship, including "plateau" behavior was reported. 

At high AK, a mild increase in fatigue crack growth rate with increasing loading frequency 

suggests that surface film destabilization promotes [H1 ABS and subsequent hydrogen ingress over 
d . Average crack tip strain rate per loading cycle (6^) is proportional to AK 7 f, where y is 
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between 2.5 and 3.5 for 7075 fatigued in moist air 004 '. As f increases at constant AK, e CT 
linearly increases and a film on the crack surface is disrupted more frequently per unit time, 
perhaps as a critical film rupture strain is exceeded (W) . If the in-tact film blocks cathodic 
hydrogen production or hydrogen atom entry to the process zone, then increased film disruption 
promotes increased hydrogen entry per unit time. That is, [H] ABS increases with increasing 
loading frequency; da/dN accordingly increases as illustrated in Figure 3a or 3b. 

The lack of a frequency effect at low AK is reasonable because e CT values are 
substantially less than those of the high AK case for the frequency range examined, and may not 
be sufficient for repeated film rupture. Surface repassivation would rapidly occur relative to the 
time interval between rupture events at the slow f CT . Hydrogen production and entry are not be 
enhanced by increasing f. For y equal to 3.0, f levels between 10 and 200 Hz are required to 
develop near-AK, h (2 MPav/ m) crack tip strain rates equivalent to those produced at higher AK 

(10 MPa'/ m) and 0. 1 < f < 5 Hz. High frequency FCP was not examined for 2090 in aqueous 
chloride. At such high frequencies, dislocation transport will not limit da/dN, however, surface 

electrochemical reactions may become important and further complicate the effect of loading 
frequency. A corollary to this argument is that increasing AK at constant frequency will strongly 
increase e CT , promote film instability, increase [H] ABS and thus increase da/dN. The effect of 
stress intensity on crack surface hydrogen concentration was not considered in the development 
of the process zone model prediction of da/dN-AK; Equation 8. 

The scenario proposed for alloy 2090 is relevant to hydrogen environmental FCP in other 
2000 series alloys that exhibit a similar frequency dependence 0 *. The situation is more complex 
for 7000 series alloys where da/dN substantially decreases with increasing frequency. In this case 
surface reaction rate and surface film stability may both govern the frequency response, with the 
former dominating surface hydrogen production. There is no evidence to suggest that hydrogen 
diffusion is growth rate limiting for 7000 or 2000 series aluminum alloys, consistent with the 
process zone model. 


V. CONCLUSIONS 

1. Environmental fatigue crack propagation rates and microscopic cracking modes, 
established in Parts I and II for Al-Li-Cu alloy 2090 in moist gases and aqueous 
chloride 0,4 *, are explained by a crack tip process zone model based on hydrogen 
embrittlement. The macroscopic growth rate (da/dN ENV ) is equated to discontinuous crack 
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advance over a distance (Aa) determined by mobile dislocation transport of hydrogen at 
plastic strains above a critical value; and to the number of load cycles (AN) required to 
sufficiently hydrogenate process zone trap sites that fracture due to a local hydrogen 
concentration-tensile stress criterion. 


The process zone model explains AK and dissolved hydrogen-dependent fatigue crack path 
transitions. For Al-Li-Cu exposed to NaCl or moist gases at low AK, FCP is along {100} 
planes due to H-induced lattice decohesion or hydride formation, and because the 
embrittlement process zone is smaller than the subgrain size; a single power-law da/dN 
versus AK relationship results. When AK increases to the level where the process zone 
just includes subgrains and dislocations reach such boundaries, environmental FCP 
transitions to the intersubgranular mode because of strong hydrogen trapping at subgrain 
boundary sites including T, precipitates; a transition occurs in the slope of da/dN-AK. A 
second transition occurs at higher AK when substantial slip plane cracking occurs due to 
mechanical fatigue damage. 


3. 


The process zone model semiquantitatively predicts the shape of the da/dN-AK relation- 
ship. For any environment, the FCP law is determined by the AK dependency and 


•oportion (©,) of each parallel microscopic cracking mode according to: 

.'i , ^ iimH ^ w 4- 


ds/dNcvio ® 


AK 2 ^ + e„*nAK (2 + 2 *-< l00}) + 0socAK (2 + 2 *- SGC) 


and including "plateau" behavior. These exponents are predicted based on the plastic 
strain range dependencies of Aa and AN for each damage mechanism. Absolute values 
of the exponents and the preexponential coetticients are not independently predicted, based 


on FCP data for alloy 2090: 

da/dN Moi „ Air = e {m) (3xl0' 9 )AK 4 ° + 0 {l()()) (lxlO R )AK 5 ' + 9 soc (2xl0 7 )AK‘ 4 
da/dN Amxlic N „ci = 0{ioo>0- 7x1 ° 7 )AK 5 ' 4 + 6 SGC (7xl0 7 )AK 2 7 
da/dN lnen = e {m} (3xl0 9 )AK 4 ° where 0 {m > = 1.0 
0 ; values, from fractographic analyses for each cracking mode, are in good agreement 
with these equations. 


The process zone model does not quantitatively predict da/dN ENV values because the local 
failure criterion for hydrogen cracking is not independently determined, the very near tip 
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plastic strain distribution is not quantitatively defined, and crack tip dislocation structures 
are not well characterized. The mechanisms by which hydrogen promotes {100} and 
intersubgranular cracking are speculative. 


5. Crack growth rates increase and the stress intensity required for the microscopic mode 
transition decreases as the environmentally produced crack surface hydrogen concentration 
increases. Crack surface films, produced by cathodic polarization in NaCl or 0 2 addition 
to water vapor, reduce hydrogen uptake and da/dN. Bulk hydrogen diffusion in the 
process zone can augment dislocation transport for slow rate FCP. 

6. For 2090 (and other 2000 series alloys) in NaCl, da/dN increases with increased loading 
frequency because increased crack tip strain rate destabilizes otherwise protective surface 
films, and enhances hydrogen uptake in the process zone. Fast crack surface chemical 
reactions, short process zone diffusion distances and rapid dislocation transport of 
hydrogen preclude rate limited cracking for frequencies less than about 150 Hz. 
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VII. APPENDICES 


A. Determination of the Coefficient a: 

The coefficient, a, is described in Equation 4 and must be determined to calculate the 
crack advance distance, Aa. It is only possible to estimate Ae crit and thus a. Assuming that a 
critical plastic strain range (Ac p ) of between 0.2% and 1% is necessary for significant crack tip 
plasticity, a will vary from 0.8 to 0.4 for peak aged alloy 2090 ( a ys = 500 MPa and E = 81.5 
GPa* 3 *). In principle Ae gri , can be determined from the implication that the threshold AK^ for FCP 
represents the condition where the peak plastic strain range is just sufficient for damage based on 
dislocation activity. AK th may be related to peak Ae p through either Equation 3, or based on an 
extrapolation of SEM measurements of crack tip strains* 

Ae p = K„ (AK - AK lh *) r (Al) 

where K,„ r and AK lh * are curve fit constants for a given alloy and environment. The singularity 
in strain as r approaches zero at the crack tip confounds the use of Equation 3, however, Ae crit 
equals 0.4% for an assumed threshold of 3.5 MPa\/ m and an r of 2 ^m (a ys = 500 MPa and E 
= 81.5 GPa). This r, the spatial resolution of the SEM method, is selected for comparison with 
the prediction of Equation Al and this threshold is a typical value for FCP in 2000 and 7000 
alloys in vacuum (3,5,45) . For 7075 in vacuum, K (> = 1.3 x 10 4 , AK^ = 0 and r = 3.0 in 
Equation Al; the corresponding plastic strain range equals 0.6% at a AK^ of 3.5 MPav/m (8 ' 78) . 
From Equation 4, a equals 0.6 and 0.5, respectively, for these two estimates of the critical plastic 

strain. 

These estimates are not rigorous. The SEM measurements for Equation Al are variable, 
do not consider crack closure effects on AK, are limited to the specimen free surface, must be 
extrapolated to yield the strain level at threshold, and are likely to depend on the fatigue cracking 
mode. The use of a vacuum AK lh in Equation 3 is justified only if it corresponds to the threshold 
dislocation activity for hydrogen accumulation. The use of AK^ for water vapor or chloride 
bearing environments is unreasonable because bulk hydrogen diffusion may be dominant for low 
crack growth rates and small process zone sizes typical of near-threshold FCP. This is the reason 
why AKtf, decreases with increasing crack surface hydrogen content, while the dislocation 
transport concept implies that AK th and Ae cril are environment independent. 
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B. Analysis of Crack Tip Opening Stress: 

Tensile stresses, <r N , normal to the crack plane and decreasing with increasing distance, 
r, from the crack tip may govern interfacial and "cleavage" cracking modes. Continuum analysis 
shows that <r N at maximum load is uniformly high (equaling 3 to 4 times cr ys ) over a distance of 
about 10 CTOD nuix l40) . Given that: 

CTOD max = r K ma > ys E (Bl) 

a yield strain of 0.6% (<7 ys /n for 2090) and T equalling 0.5 6 , the distance for high o N can be 
restated in terms of d cz : 


10 CTOD inax = {1.2/Ml - R) 2 ]} d ez (B2) 

When a equals 0.5, <r N is high over 3 d c/ at a stress ratio of 0. 1 and over 250 d ez for R equaling 
0.9. This latter value is typical of the FCP experiments conducted with alloy 2090 (3) . 

C. Determination of the Relationship Between AN and AK: 

In the process zone model, crack advance is presumed to occur over a distance, Aa, when 
AN equals the number of load cycles sufficient to achieve the critical trapped hydrogen 
concentration within the crack tip process zone, d tf/ . To determine the AK dependence of AN, 
assume that AN inversely depends on the total number of mobile dislocations within the 
embrittlement zone volume, N m , raised to some power, -<f>: 

AN = cN./ (Cl) 

where f is a constant. Since mobile dislocation density (p) is proportional to plastic strain* 96 *, N m 
equals the product of p, d ez and unit process zone thickness. For the more complex case, where 
A( p depends on distance ahead of the crack tip, N m is described by the integral of the plastic 
strain range distribution from the crack tip to d c/ , per unit process zone thickness. That is: 


6 The proportionality constant, F, varies between 0.4 (finite element) and 1.0 (Dugdale strip yield) depending on the 
method of crack tip field analysis* 40 *. 
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(C2) 


« d cpz 

N m = ^ Ae p (x) dx 

The plastic strain distribution ahead of the crack tip is given in Equation 3. x is a small distance 
from the crack tip were the continuum analysis predicts a strain singularity, but within which the 
actual strain is finite and perhaps a constant. Substituting the expression for Ae p (x) into Equation 
C2 and integrating yields: 

N m = [w<r y yE][{(n+l)a W(n+l) d cp/ /n}{l-(na ,/(n+,) /(n 4- 1))} - 

{{(n+ |)d cpz l/(n+1) x*" /(n+1> }/n} + x*] (C3) 

For low work hardening materials such as Al-Li-Cu, n is on the order of 0.1, for a = 0.5 from 
Equation 4, and for o> = 1.0‘, Equation C3 is simplified to. 

N m - [a y J E] [0.8d cp/ - 1 1 .Ox* 0,w d cp/ 09 + x* (C4) 

Since d cp7 is proportional to AK 2 (Footnote 2), and tor small x on the order of 0.05 times d cpz . 

N m = 0.02AK 2 /<r ys E (C5) 

Combining Equations C5 and Cl yields the AK dependence of AN for a low work hardening 
material (n =0): 

AN = g'a ys *E*/AK 2 * (C6) 
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Cycles required for crack extension (AN) 


Alloy 7075T651 



Stress intensity range (MPaVln) 


Figure 1 A comparison of calculated load cycles for discontinuous crack growth in alloy 
7075 exposed to either moist air or inert environments. AN values are based on 
either the process zone model with measured da/dN, or SEM measurements of 
crack tip strain range with a Coffin-Manson failure criterion <8) . 
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I'igure 2 A schematic drawn to scale showing the crack tip process zone at constant AK for 
hydrogenous environments; (a) Regions B and C SGC where the embrittlement 
process zone is larger than subgrain size, and (b) low AK {100} FCI\ Region A 
in Fable 1, where the process zone is smaller than the subgrain size. 
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Log da/dN Lo 9 da/dN 



Log AK 


Figure 3 Da/dN versus AK relationships based on the process zone model for hydrogen 
enhanced FC1 > in Al-Li-Cu alloys, (a) and (b) For a transition between single 
inodes of cracking, [H] ABS affects FCP for dislocation transport of hydrogen within 
d c7 and for dislocation transport augmented by bulk hydrogen diffusion, (c) 
Complex da/dN-AK relationships result when multiple FCP modes operate in 
parallel at fixed [H] ABS . 
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Fatigue crack growth rate (mm/cycle) 


Alloy 2090 
Moist air, 5 Hz 
Kmax = 17 MPa VTn 


L-S orientation 
L-T 

L-T, constant AK, R = 0.54 
L-S, constant AK, R = 0.60 
da/dN = 3 x 10 ^ 1 1 } AK^ - ^ + 10 ^ ^{1 00 

2 x 10' 7 Ooz-i- AK 2 * 4 
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Stress intensity range (MPa V~m ) 


A comparison of measured and predicted da/dN-AK ''plateau" behavior for al!o\ 
2090 in moist air. 
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Figure 6 
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A schematic depicting: (a) occluded cell electrochemical reactions in Al-Li-Cu 
leading to crack surface hydrogen production, [H] Ans , and (b) crack tip pH and 
atomic hydrogen production rate as a function of electrode potential. 
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